Mon. Not. R. Astron. Soc. 000, [THg] (2008) 



Printed 11 February 2010 



(MN M?eX style file v2.2) 



On the Nonthermal Emission from the Supernova 
Remnant W51C 



o 

(N 
(X, 



6 

o3 



> 

(N 
(N 

(N 
O 

o 



Jun Fang , Li Zhang * 

Department of Physics, Yunnan University, Kunming, China 



11 February 2010 



ABSTRACT 

The middle-aged supernova remnant (SNR) W51C is an interesting source for the 
interaction of the shell with a molecular cloud. The shell emits intense radio syn- 
chrotron photons, and high-energy gamma-rays from the remnant have been detected 
using the Fermi Large Area Telescope (LAT), the H.E.S.S. telescope, and the Milagro 
gamma-ray observatory. Based on a semi-analytical approach to the nonlinear shock 
acceleration process, we investigate the multiband nonthermal emission from W51C. 
The result shows that the radio emission from the remnant can be explained as syn- 
chrotron radiation of the electrons accelerated by a part of the shock flowing into the 
ambient medium. On the other hand, the high-energy gamma-rays detected by the 
Fermi LAT are mainly produced via proton-proton collisions of the high-energy pro- 
tons with the ambient matter in the molecular cloud overtaken by the other part of 
the shock. We propose a possible explanation of the multiband nonthermal emission 
from W51C, and it can be concluded that a molecular cloud overtaken by a shock 
wave can be an important emitter in GeV 7-rays. 
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1 INTRODUCTION 

SNRs are broadly thought as the primary accelerators of 
the Galactic cosmic rays accelerated th rough diffusive shock 
acceleration in the shock waves (e.g., iBIandford fc Eichlerl 
ll987l : lBerezhko fc Volkll2007l ). Observations of SNRs in the 
high-energy 7-ray band provide evidence for acceleration 
of par ticles to relativistic energies (e.g., lAharonian et all 
l2008ch . However, both electrons and nuclei can be accel- 
erated to relativistic energies, and the high-energy 7-rays 
from SNRs can usually be explained either as proton-proton 
(p-p) interaction between the relativistic protons and the 
ambient matter or as inverse Com pton scattering of the 
electrons on the soft photons (e.g., Yamazaki et al.l 20061: 
Aharonian et all 120071; iBerezhko fc Volkl |200S|; iFang et ail 



20081 . 120091 : IZirakashvili fc Aharon ian 2010). The interpre- 



tation of the 7-rays produced via the inverse Compton scat- 
tering makes the evidence of the argument that the nuclei in 
the cosmic rays are accelerated by the SNR shells indirectly. 

P-p collisions can be greatly enhanced if a SNR shell 
interacts with a molecular cloud, and observations of 7-rays 
from this system provide evidence of the nuclei in the cos- 
mic r ays being accelerated by SNRs (e.g.. [Aharonian et al.l 
1 19941 ). Moreover, 7-rays from such systems have probably 
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been detected with H.E.S.S. (e.g., lAharonian et ai]|2008al lbl) 
and EGRET (|Hartman et all 1 19991 '). 

The SNR W51C appears as a partial shell of ~ 30' 
diameter in the radi o continuum with an open northern part 
(|Moon fc Kodl 19941 ) . From the observations of the two 1720 
MHz OH masers toward the north of the shell (|Green et al.l 
Il997t) and the detection of shocked atomic and molecular 
gases in this direction (|Koo fc MoorJll997al lbh. it is argued 
that this remnant is interacting with a molecular cloud. 

Both shell-type and center-filled morphologies were in- 
dicated in the X-ray observations with ROSAT, and an age 
of ~ 3 x 10 4 yr and an explosion energy ~ 3.6 x 10 51 erg 
were estimat ed base on eithe r the Sedov or the evapora- 
tive models (|Koo et al.l [T995) . An extended source, CXO 
J192318. 5+1403035, composed of a relatively bright core 
surrounded by a diffuse envelope was discovered with ASCA 
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20021 ). and it was also detected with Chandra 



2005). The core-envelope structure and its spec 



tral properties suggest that the source is a pulsar wind neb- 
ula (PWN) associated with the SNR 51C, and a putative 
pulsar is argued to exist inside it. 

In TeV 7-rays, an extended source HESS J1923+141 
coincident with W5 1C have bee n found us i ng H.E .S.S. 
(|Fiasson et al.1 l2009h . Moreover, lAbdo et al] ( |2009ah re- 
cently reported a possible excess of multi-TeV gamma-rays 
towards W51C with a flux of 39.4 ± 11.5 x 10~ 17 TeV" 1 
cm -2 s _1 at 35 TeV from the Milagro observation. Further- 
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more, the GeV 7-rays from the SNR, which is known to 
be interacting with a molecular cloud, with a luminosity of 
~ 1 X 10 36 erg s _1 for a distance of 6 kpc, have been de- 
tected wit h the LAT on board the Fermi Gamma-ray Space 
Telescope (|Abdo et al.ll2009bh . The 7-ray spectrum cannot 
be fit ted with a single po wer law and steepens above a few 
GeV. [Abd o ct al. (2009b) discussed the nonthermal radia- 
tive properties of the radio emission and 7-rays from the 
SNR using a broken power-law with the same index for the 
momentum distribution of the radiating electrons/protons. 
Their result shows that p-p interaction is more likely the 
main process to produce the observed 7-rays since the ob- 
served radio synchrotron spectrum cannot be well repro- 
duced in the bremsstrahlung-dominated case and an origin 
of inverse Compton scattering for the 7-rays needs a low 
density of hydrogen (< 0.1 cm -1 ) in the SNR, which con- 
flicts with the constraints from observations in the radio and 
X-ray bands. 

Semi-analytical methods to the nonlinear diffusive 
shock acceleration process have been widely used to inves- 
tigate the nonthermal radiative properties of SNRs (e.j 
iMorlino et all 120091 ; iFang et all 120091 ; lEllison et ail I201C 
In this paper, we study the nonthermal radiative proper- 
ties of the SNR W51C based on a semi-analytical method 
to the nonlinear diffusive shock acceleration mechanism. In 
our model, a part of the SNR shell transports into the molec- 
ular cloud (MC), then the shock is greatly decelerated due 
to the large density of matter in the MC. As a result, the 
spectra of the accelerated protons and electrons are deter- 
mined in a case with a low Mach number. We find that the 
observed 7-ray spectrum for W51C can be well reproduced 
via p-p collisions of the accelerated protons in the part of 
shell interacting with the MC, and the radio emission from 
the SNR can be explained as the radiation of the electrons 
accelerated by the other part of the SNR shock freely ex- 
panding into the ambient interstellar space. 



2 MODEL AND RESULTS 

The pitch-angle averaged steady-state distribution of the 
protons accelerated at a shock in one dimension satisfies 
the diffusive transport equation (e.g.. | Malkov fe Drury|200ll ; 
lBIasill2002l ; lAmato, Blasi fe Gabicill20Qgft . 



d_ 

d.v 
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df(x,p) 
dx 

1 du df(x,p) 
3 dx^ dp 



+ Q(x,p) = 0, (1) 



where the coordinate x is directed along the shock normal 
from downstream to upstream, D is the diffusion coefficient 
and u is the fluid velocity in the shock frame, which equals 
M2 downstream (x < 0) and changes continuously upstream, 
from ui immediately upstream (x = + ) of the subshock to 
uo at upstream infinity (x — +00). For the Bohm diffusion, 
D — pc 2 1 (3eB) , where B is the local magnetic field strength. 
With the assumption that the particles are injected at im- 
mediate upstream of the subshock, the source function can 
be written as Q(x,p) = Qo(p)S(x). For monoenergetic injec- 
tion, Qo(p) is 



Qo(p) 



T) n gaSi imi 



<*(P-Pinj) 



(2) 



where pi n j is the injection momentum, n gaSl i is the gas 
density at x = + and r\ is the fraction of particles 
injected in the acceleration process. With the injection 
recipe known as thermal leakage, 77 can be described as 
r? = 4(R Bub - l)g 3 e~ s2 /37r V2 dBlasi. Gabici fe Vannoiiil 



120051 ; lAmato. Blasi fc Gabicil 120081 ), where ik„ h = u%/u 2 
is the compression factor at the subshock and £ is a pa- 
rameter of the order of 2-4 describing the injection mo- 
mentum of the thermal particles in the downstream region 
(Pmj = £Pth,2). We use £ = 3.5 as in lAmato. Blasi fc Gabicil 
(2008), Pth,2 — (2m p fcBT2) 1//2 is the thermal peak momen- 



tum of the particles in the downstream fluid with tem- 
perature T2, m p is the proton mass and fee is the Boltz- 
mann constant. Assuming the heating of the gas upstream 
is adiabatic, with the conservation condition of momentum 
fluxes between the two sides of the subshock, we can de- 
rive the relation between the temperature of the gas far 
upstream T and T 2 , i.e., T 2 = (•y g Mo/Rtot)[Rsu.b/Rtot - 
1/iitot + (l/7 9 A/,?)(^tot/J?sub) 79 ]r , where M is the fluid 
Mach number far upstream, 7? to t = Uo/uz is the total com- 
pression factor, 7 S is the ratio of specific heats (7 9 = 5/3 for 
an ideal gas) 

Particles with larger momenta move farther away from 
the shock than those with lower momenta, only particles 
with mome ntum ^ p can reach the point x p (see details in 
lBlasill2002D . and thus the pressure of the accelerated parti- 
cles at the point can be described as 

Pcr.p = y ^ Pmax dp' P ' 3 v(p')Mp') , (3) 

where v(p) is the velocity of particles with momentum p. 
Furthermore, the particle distribut ion function fo(p) at the 
shock can be implicitly written as (|Blasil |2002| ) 



Mp) 



3.Rti 



RtotU(p) 



x exp 



dp' 3R tot U(p') 



p> Rt ot U(p>) - 1 



(4) 



where ni sm is the gas density far upstream (x = +00), 
Pmax is the maximum momentu m of the accelerated pro- 
tons, and it can be estima ted as ijLagaee fc Cesarskvlll983l : 
lAharonian fc Atovanlll999l ) 



B 



Pm - = 5 " 3X 104 iTO^H 108cm s-0 mpC ' (5) 

for a "Bohm limit" according to the standard diffusive shock 
acceleration theory, where B the magnetic field strength, ts 
is the age of the host SNR in units of 10 3 yr. U(p) can be 
solved using an equation deduced from the conservation of 
the mass and momentum fluxes with the boundary condition 
U(pinj) = -R S ub/#tot and (7(p max ) = 1, and then a value of 
-R su b can be achieved by an iterat ive procedure to satisfy the 
boundary conditions (|Blasill2002l ). 

The electrons have the same spectrum of the protons 
up to a maximum energy determined by s ynchrotron losses , 
and the spectrum can be described as (e.g.. lFang et al.ll2009l ) 



f e (x,p) = K ep f(x,p)exp(-E(p)/E m ^e), 



(6) 



where E(p) is the kinetic energy of the electrons, and 
the electron/proton ratio K ep is treated as a parameter. 
By equating the synchrotron loss time with the acceler- 
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ation time, the cutof f energy E n 
jBerezhko et alj|2002f ) 



can be estimated as 



E n 



6x10' 



"o 



10 8 cm s _1 



-Rtot(l + ^?sub-Rtot) V ^ 



Assuming the acceierated particles distribute homoge- 
neously and most of the emission is from downstream of 
the shock, and using the distribution function at the shock 
to represent the particle distribution in the whole emitting 
zone, the volume-averaged emissivity for photons produced 
via p-p interaction can be written as 



Q{E) = 47rn gas / dE p J p {E p 



da{E, E p 
dE 



(8) 



where E p is the proton kinetic energy, n gas is the ambi- 
ent gas number density, and J P (E P ) = vp 2 fo(p)dp/dE p is 
the volume- averaged proton density and v is the particles' 
velocity. We use the diffe rential cross-s e ction for photons 
da(E, E p )/dE presented in lKamae et al.l (|2006h to calculate 
the hadronic 7-rays produced via p-p collisions. Finally, the 
photon flux observed at the earth can be obtained with 



F(E) 



VQ(E) 
4%d 2 



(9) 



where d is the distance from the earth to the source and V 
is the average emitting volume of the source and can be es- 
timated by V « (47r/3)-Rf nr /-Rtot l|Ellison et alJl200Ch . here 
i? S nr is the radius of the SNR. 

The SNR W51C has an age of ~ 3 x 10 4 yr |Koo et all 
120051) and a radius of ~ 30 pc for a distance of 6 kpc 
jMoon fc Kool 1 19941 ). With d = 6 kpc, i? snr = 30 pc, 



ton spectra and the spectral energy distribution of the p-p 
collisions for the protons with different Mach number are 
illustrated in Fig[T] The results are weakly sensitive to the 
temperature upstream of the shock, and a temperature of 
0.5 x 10 6 K is used in the calculation, which corresponds 
to a sound spe ed of 83 km s -1 . T he observe d results with 
the F ermi LAT (lAbdo et al.ll2009bl), H.E .S.S. (|Fiasson et all 
l2009f ). and Milagro (|Abdo et all |2009ah are shown in the 
figure for comparison. For each Mo, rii sm is normalized to 
make the resulting flux at 9 GeV consistent with the one 
observed with the Fermi LAT. We find that the spectrum 
given by the Fermi LAT can only be reproduced with a low 
Mach number ~ 3 with those parameters. For such a low 
Mach number, the reaction of the accelerated particles on 
the shock is negligible, and the value of i? su b deviates lit- 
tle from that of i?tot- In fact, for Mo = 2.8, -R su b and i?tot 
are 2.83 and 2.91, respectively. The index for the momen- 
tum distribution of the accelerated protons/electrons can 
be estimated as 3(7 S + l)Mrj/(2(Mrj — 1)), which is equal 
to 4.5 for Mo = 3 and 7 3 = 5/3. With such an index, the 
observed radio emission cannot be well explained with the 
synchrotron emission of the electrons, although the emission 
observed with the Fermi LAT in the 7-ray band can be well 
reproduced as the p-p interaction of the accelerated protons . 

Observations of the OH masers (|Green et at! Il997l) 
and th e shocked atomic and molecular gas ( Koo fc Moonl 



Il997al lbl) in the direction of W51C indicate that a part of 
the SNR shell is interacting with the MC. As a result, the 
velocity of this part of the shell decreases significantly, and 
the observed 7-rays can be reproduced as p-p collisions ac- 
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Figure 1. The resulting proton spectra (upper panel) and the 
spectral energy distribution of photons from the p-p collisions 
(lower panel) for Mq = 15, n; sm = 0.34 cm~ 3 (solid line); Mq = 
10, n; sm = 0.53 cm -3 (dashed line); Mq = 5, n; sm = 1.51 cm -3 
(dotted line); and Mo = 2.8, n; sm = 14.4 cm 
line). The other parameters are d = 6 kpc, i? S nr 
0.5 xlO 6 K, and Pmax — 10^ m p c. For each Mo, Tii sm is normalized 
to make the resulting flux at 9 GeV consist with the one observed 
by the Fermi LAT. The flux points for the SNR W51C obtained 



3 (dash-dotted 
= 30 pc, T up = 



with the Fermi LAT llAbdo et al. 



H.E.S.S. l lFiasson et al.l 
2009aJ) are shown in the figure 



120091) . and Milagro l lAbdo et alj 
for comparison. 



celerated in this part of the shell due to its low Mach number. 
The other part of the shell can freely expand into the inter- 
stellar medium with a relatively low density. Therefore, we 
assume about / = 0.4 of the SNR shell transports into the 
MC with a Mach number of 2.8 now, and the other part of 
the shell have a relatively large Mach number, Mo = 10, cor- 
responding to a shock speed of 830 kins -1 for T up — 0.5 xlO 6 
K. The synchrotron photons in the lower energy band are 
attenuated via free-free absorption in the journey from the 
source to the earth, and the absorption co efficient can be 
described as (|Ginzburg fc Svrovatskiill 19651 ) 



as H = 10" 



v 2 T t 



3/2 



17.7 + In 



T, 



3/2 



(10) 



where v is the electron frequency, n e and T e are the density 
and temperature of the intervening ISM, respectively. We 
use T e = 8000 K and n e d = 9.3 x 10 22 cm -2 for W51C. The 
resulting multiband nonthermal emission for the two parts 
of the shell are shown in Fig[2] The soft photons involved 
in the inverse Compton scattering include the cosmic mi- 
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Figure 2. The resulting synchrotron radiation (dashed line), 
bremsstrahlung (dotted line), inverse Compton scattering (dash- 
dotted line) and p-p interaction (solid line) for the shell expand- 
ing into the low-density interstellar medium (upper panel) and for 
the shell interacting with the MC (lower p anel) with / = 0.4. The 
observations in radio jMoon fc Kodll994f) and high -energy 7-rays 
llAbdo et al., 2009a b; Fiasson ct al. 2 009T ). The other parameters 
for the upper panel are d = 6 kpc, i? snr = 30 pc, T up = 0.5 X 10 6 



K, M = 10, 



0.08 cm" 



2 x 10" 



13 



150 



/iG. Those for the lower panel are Mo = 2.8, 

^ism — 21.4 cm , 

B = 15 fiG, and the others arc the same as the upper panel. In 
the lower panel, rii sm is determined by making the resulting flux 
at 9 GeV consist with the one observed by the Fermi LAT, and 
the error of n; sm is n; sm = 6.0 cm -3 due to the uncertainty in 
the observation. 



crowave background, infrared background light and optical 
star light, with the energy densities for the three component 
0. 25, 0.9 and 1 .0 MeV cm -3 , respectively, similar to those 
in lAbdo et all (|2009bl ). 

The observed emission with a shell-type morphology in 
the radio band from W51C can be explained as the syn- 
chrotron radiation from the electrons in the shell not en- 
countered by the MC with a Mach number of 10. A mag- 
netic fi«ld_st : ren|[Ui^^00jiG for SNRs 
(e.g JParizot et alfeOOfjlTanaka et al.ll2008l ). With B = 150 
/iG, the maximum energy of the protons is ~ 1.6 x 10 3 
TeV (Eq[5]) for an age of 3 x 10 4 yr, which is near the 
"knee" energy of the Galactic cosmic rays. With a den- 
sity of ni sm = 0.08 cm -3 , the resulting flux from p-p in- 
teractions at 35 TeV is consistent with the observation with 
Milagro. The electron to proton ratio K ep cannot be the- 



oretically determined now, which is usually treated as a 
parameters limited by the comparison of the resulting flux 
from the model with the observations for a given source 



(e.g., iBerezhko et al.ll2002l; rMorlino et all l2009l ; iFang et all 
|200S| . 120091 ; IZirakashvili fc Aharonianl I2010T ). Its value can 



be from ~ 10 4 for young SN Rs to ~ 10 ^ for old ones 
|Zirakashvili fc Aharonianl I2010T ). In our model with Mo = 
10 and B = 150 ^G for the part of the shell of W51C not 
colliding with the MC, the radio must be 2 x 10~ 2 to have 
the detected radio emission explained as the synchrotron ra- 
diation of the electrons, which is about twice the value in 
the galactic cosmic rays at 10 GeV around the earth. In or- 
der to produce the observed fluxes in the radio bands, the 
ratio K ep cannot be very smaller than 0.02 except a much 
stronger magnetic field is used. In such a case, the emis- 
sion above 1 TeV is dominated by p-p interactions of the 
accelerated protons containing a kinetic energy of 3.4 x 10 50 
erg, and the TeV photons at 35 TeV observed with Mila- 
gro are from this mechanism (see the upper panel in Fig(2|. 
The maximum energy of the accelerated electrons is 1.7 TeV 
(Eq[7]). With the free- free absorption taken into account, the 
detected fluxes in the radio band can be reproduced with the 
synchrotron radiation of the accelerated electrons. 

For the part of the SNR shell interacting with the MC, 
with Mo = 2.8, rii sm = 21.4 cm -3 is needed to make the 
resulting flux of p-p collisions consist with that at 9 GeV de- 
tected with the Fermi LAT. The flux from the accelerated 
electrons is sensitive to the electron to proton ratio K ep , 
which is set to 2 x 10~ 2 as for the part not interacting with 
the MC. With this ratio, the magnetic field downstream of 
the shock cannot be much larger than 15 fiG to avoid the 
resulting synchrotron flux to conflict with the observed flux 
points in the radio band. In such a case, the maximum ener- 
gies of the electrons (15 maXie ) and the protons (£ max ) are 3.0 
and 12 TeV, respectively. With these reasonable parameters, 
the high-energy 7-ray spectrum of the SNR W51C observed 
with the Fermi LAT and H.E.S.S. can be well reproduced 
(see the lower panel in Fig[2]). These 7-rays are mainly pro- 
duced via p-p interaction between the protons in the shell 
interacting with the MC and the ambient matter, and the 
kinetic energy contained in the protons is about 5 x 10 50 erg. 

There are uncertainties on the parameters we used to 
reproduce the multiband observed flux points for the SNR 
W51C in our model. For instance, the value of the ambi- 
ent density is sensitive to the Mach number and the whole 
number of the accelerated particles, which is determined by 
the emitting volume with a known spectra of the acceler- 
ated protons. In this paper, the ratio of the shell colliding 
with the ambient MC is assumed to be / = 0.4 based on 
the observed extension of the 7-rays is similar as that in the 
radio continuum. For the part of the shell interacting with 
the remnant, a Mach number of 2.8 is derived to make the 
resulting fluxes from p-p collisions consist with the shape of 
the flux obtained with the Fermi LAT. With such a Mach 
number, a density of ni sm = 21.4 cm -3 is needed to repro- 
duced the observed flux in the 7-ray band. The distribution 
of the protons is oc ni sm , and this relation is also valid for 
the rate of p-p interaction. As a result, we can only limit 
the value /n 2 sm from equaling the resulting flux of p-p in- 
teractions to the observed one, and the density varies from 
25 cm" 3 for / = 0.3 to 19 cm" 3 for / = 0.5. Moreover, the 
error of the observed 7-ray fluxes is about 50%, and then the 
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uncertainty of n- m 
for / = 0.4. 



limited in our model is about 6.0 cm 



3 DISCUSSIONS AND SUMMARY 

W51C is an interesting source for the interaction of the 
shell with the MC from the observations of OH masers and 
shocked atomic and molecular gases in the direction of the 
source. It is also the first SNR detected with the Fermi LAT 
in the GeV band. The SNR appears a shell-type morphology, 
whereas the GeV signals are clumpy and extended within the 
radio shell. The 7-rays steep above several GeV, and the 
spectrum cannot be reproduced via p-p collisions with an 
index of ~ 2.0 for the relativistic protons. Cosmic ray elec- 
trons with a power -law spectrum ca n be easily steepened by 
unity duo to losses (|Kardashevlll962l 'l. and it is plausible that 
this steepening can be explained in the scenario of a power- 
law electrons encountering significant losses. However, the 
multiband nontherma l radiative p r opertie s of W51C have 
also been studied by lAbdo et al.1 (|2009tJ ) using a broken 
power-law spectrum for the energy distribution of the elec- 
trons/protons. They found the observed flux points in the 
radio and GeV bands cannot be well reproduced by assum- 
ing the radio emission and high-energy 7-rays are produced 
by the same population of electrons, and only in the pion- 
decay scenario in which the 7-rays are produced by the pro- 
tons can the multiband spectrum be well fitted. Moreover, 
it seems unlikely that the GeV s ignals are from t he PWN 
powered by the putative pulsar (|Koo et alj|2005h because 
the nebula is much smaller with a few pc extension in X- 
rays compared with the GeV sourc e and its radio emission 
is much weaker l|Abdo et al. I l2009bl ). 

We model the source in the case that one part of the 
shell transports into the MC with dense matter and the 
other part expands into the interstellar space with a rel- 
atively small density. The spectra of the particles in the 
two parts of the shell are treated in a semi-analytical ap- 
proach to the nonlinear shock acceleration. Note that in this 
paper the shock structure is assumed to be quasi-parallel 
rather than quasi-perpendicular. The maximum energy of 
the protons accelerated in the quasi-perpendicular shocks 
can be significantly higher than that in the quasi-parallel 
cases with a much smaller perpendicular diffusion efficient 
compared with the parallel one, and a q uasi-parallel shocks 
cannot possibly be true over all of An |jokipiil |l987t). The 
contribution of the secondary e pairs produced from p-p 
interactions on the nonthermal emmission is not taken into 
account in this paper. Including it needs a time-dependent 
model to describe the evolution history of the SNR and the 
collision between the part of the shock and the MC, and 
this complicated process could bring more assumptions in 
the model; moreover, with K ep — 2 x 10~ 2 , the emission of 
the secondary pairs is usually insignificant compared with 
that from the accelerated electrons. In fact, the multiband 
observed spectrum can be reproduced with our model, thus 
we do not take into account the emission from the secondary 
e pairs produced in p-p collisions. Our results indicate that 
the 7-rays detected with the Fermi LAT and H.E.S.S. are 
mainly from the p-p collisions of the protons in the shell 
interacting with the MC, whereas the radio emission are 
produced via synchrotron radiation of the electrons accel- 



erated in the shell not interacting with the MC; moreover, 
the high-energy 7-rays obtained with Milagro at 35 TeV are 
possibly from p-p interactions of the protons accelerated in 
the shell not interacting with the MC. 

A MC can be illuminated by the particles escaped from 
a nearby SNR, and high-energy 7-rays can be produced 
when these particles d iffuse into the cloud |Gabici et al.l 
120091 ; iFuiita et al.l [20091 '). For example. Pforres et al.1 (|2008l ) 
interpreted the source MAGIC J0616+225 as TeV emission 
from the particles escaped from the SNR IC 443 and diffus- 
ing into a MC located about 20 pc from the SNR. Moreover, 
the nonthermal radiative properties of the the open cluster 
Westerlund 2 and the old-age mixed-morphology SNR W28 
can also be modeled as MCs illuminated by hi gh-energy 
particles escaped from SNRs (|Fuiita et al.l 120091 ). On the 
other hand, 7-rays can also be effectivel y produced from 
a shock wave colliding wi th a MC (e.g., lAharonian et al] 
1 19941 ; IZhang fe Fanell2008h . 7-ray emission produced from 
a MC overtaken by a shock wave with a low Mach num- 
ber has also been used to discuss the emission properties 
of t he TeV sources HES S J1745-303 (A) and HESS J1714- 
385 (|Fang fe Zhandl2008l '). In this paper, we model the SNR 
W51C as the MC around the remnant overtaken by part of 
the shell, and the other part of the shell expands into the rel- 
atively tenuous interstellar medium. The spectrum obtained 
with the Fermi LAT can be reproduced as 7-rays from p-p 
collisions of the protons in the shell with a low Mach num- 
ber. 
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